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ABSTRACT 

Recent cosmological modeling efforts have shown that a local underdensity on scales of a few hundred 
Mpc (out to z ~ 0.1), could produce the apparent acceleration of the expansion of the universe 
observed via type la supernovae. Several studies of galaxy counts in the near- infrared (NIR) have found 
that the local universe appears under-dense by ~ 25 — 50% compared with regions a few hundred Mpc 
distant. Galaxy counts at low redshifts sample primarily L ^ L* galaxies. Thus, if the local universe 
is under-dense, then the normalization of the NIR galaxy luminosity function (LF) at z > 0.1 should 
be higher than that measured for z < 0.1. Here we present a highly complete (> 90%) spectroscopic 
sample of 1436 galaxies selected in the i/— band (1.6^m) to study the normalization of the NIR LF at 
0.1 < z < 0.3 and address the question of whether or not we reside in a large local underdensity. Our 
survey sample consists of all galaxies brighter than 18*'' magnitude in the if— band drawn from six 
widely separated fields at high Galactic latitudes, which cover a total of ^ 2 deg^ on the sky. Wc find 
that for the combination of our six fields, the product <j)*L* at 0.1 < z < 0.3 is ^ 30% higher than 
that measured at lower redshifts. While our statistical errors in this measurement are on the ^ 10% 
level, we find the systematics due to cosmic variance may be larger still. We investigate the effects 
of cosmic variance on our measurement using the COSMOS cone mock c atalogs from the Millennium 
simulation and recent empirical estimates of cosmic variance derived bv iDriver fc Robothaiiil (fMol) . 
We find that our survey is subject to systematic uncertainties due to cosmic variance at the 15% level 
(1 cr), representing an improvement by a factor of ^ 2 over previous studies in this redshift range. We 
conclude that observations cannot yet rule out the possibility that the local universe is under-dense 
at z < 0.1. The fields studied in this work have a large amount of publicly available ancillary data 
and we make available the images and catalogs used here. 

Subject headings: cosmology: observations — galaxies: fundamental parameters 



1. INTRODUCTION 

The size of large-scale structures in the local universe 
(e.g. sheets, voids, and supcrclusters of galaxies), and 
our location among them, is of critical importance to the 
interpretation of observational results. Cosmic variance 
due to large-scale structure can lead to systematic vari- 
ations in observational data. Such systematics can dom- 
inate over other sources of error if a volume sufficient 
to average over cosmic variance has not been sampled. 
However, it remains unclear just what the upper limit on 
the size of large-scale structure is, and hence what volume 
constitutes a representitive sample of the universe. If the 
typical size of local large-scale structure is much greater 
than 100 Mpc, then local measurements of the Hubble 
constant, and other cosmological observables, could still 
harbor large systematic errors. 

Cold dark matter simulations that include a cosmolog- 
ical c onstant (ACDM mo dels), such as the Millennium 
Run ([Springel et al.l [20051 ) . predict that the largest dark 
matter structures in the universe should be 100 Mpc 
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in extent. Observed luminous matter large-scale struc- 
tures, most notab ly the > 400 Mpc Sloan Great Wall 
(|Gott et al.ll2005D . demonstrate the existence of struc- 
ture on larger scales than simulations predict (assum- 
ing luminous matter traces dark matter). It has also 
been shown that voids on scales similar to that of the 
Sloan Great Wall may explain the "cold spots" in the 
cosmic microwave background ()Inoue &: Silkll2006l ). It is 
not yet clear whether such structures represent extreme 
non-linearities in the matter distribution of the universe, 
or if inhomogeneities on several hundred Mpc scales are 
typical. 

Cosmologists have proposed so called "void models" 
as alternatives to ACDM, that invoke a large local 
underdensity to explain the apparent acceleration 
of the expansion of the unive rse (Alncs et al. 2006; 
Chung & Romano '2006'; 'Enavist & Mattsson 2007; 
Garcia-BeUi do & Haugb0lle 200& Al exander et ah 20091: 

[2001 iFebruarv et al l 
[Biswas et al.l [20101: 
Clarkson fc Maartend 



Garcfa-Bellido fc Haugb0lle[ 
20101: ICelerier et al.l I2OI 
Marra fc PaakkonenI [201 



20101 : [Boleiko fc Sussmanl[2011| ). The scale of the voids 



proposed in these models ranges from a few hundred 
Mpc to several Gpc in radius. The basic premise of 
void models is that if we, as observers, live near the 
center of a large underdensity, then we would witness 
a local expansion of the universe that is faster than 
the global expansion, simply due to the void being 
evacuated toward higher-density surroundings via grav- 
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ity. This would provide for a focally measured Hubble 
constant that is higher than the global value and look 
observationally like an accelerating expansion. 

In their current form, void models with A = 
appear not to be viable alternatives to ACDM, as 
they have trouble fi tting the e ntire range of cosmo- 
logical observabl e s (iGarcfa-Bell ido & Haugb0llc 2008; 
Zibin et all [200l iMoss et all 120 lit .Zhang fc Stebbind 



201 It iRiess et al.lfeoi 1[ ) . However, the exploration of this 



class of cosmological models, and other inhomogeneous 
models, has highlighted the need for a more comprehen- 
sive understanding of large-scale s tructure in the universe 
and our location within it (e.g. iMarra fc Notaril[20nl: 
IBoleiko et ail [2011 ^B^k Cliftonll20ia). In particular 
so ca l led "minimal void" scenari os (e.g.. [Alexander et al.l 
120091 iBolejko fc SussmanI 120 111 ) have shown that very 
simple models that place the observer near the center of 
a void that is ~ 250 h^^ Mpc in radius (to z ~ 0.1) and 
^ 50% under-dense compared to its surroundings are suf- 
ficient to explain the apparent acceleration observed via 
type la supernovae. While these models are simplistic, 
they point out that our location within local structure 
may have profound implications for our measurement of 
cosmological observables. 

The existence of a large local void would be consistent 
with studies of near-infrared (NIR) galaxy counts that 
have found the local space density of galaxies may be 
low by 25 — 50% compared to th e density at distances 
of ^ 300 h-^ Mpc or z ^ 1 fKeenan et al." '2010a; 
BussweU et al.ll2004HFrith et al.ir2003, ,2005^ .Huang et al.. 



1997[ ). Galaxv counts in the NIR at low reds hifts are pri- 
maril y sampling L ~ L* or brighter galaxies ()Barro et al.l 
|2009() . Thus, if the universe at z < 0.1 is under-dense, 
the normalization of the NIR luminosity function (LF) at 
z > 0.1 should be higher than locally measured values. 
In this paper we explore the normalization of the NIR 
LF as a function of redshift to test for the existence of a 
large local underdensity. 

At z < 0.1, the NIR LF has been studied using 
samples selected from large NIR surve ys, such as the 
Two Micron AU Sky Survey (2MASS. tskrutskie et aD 
\2()M\ and the UKIRT Infrar ed Deep Sky Survey 
(UKIDSS, [Lawrence et al.l [2007D . combined with spec- 
troscopy fro m redshift surveys, such as the CfA2 red - 
shift survey (IGeller fc Huchraill98llHuchra et al.llT99l . 
the Two-degree Field Galaxy Redshift Survey (2dFGRS, 
iColless et al. 2001) and the Sloan Digital Sky burvey 
(SDSS. IYork et al.ll2OOO0 . 

Using these data, the (z) ^ 0.05 NIR LF has been rela- 
tively well established (Cole et al.ll2001l : iKochanek et al.l 



ly A 

I200l][ iJones et al.ll2006D . Studies at slightly higher me- 
dian redshifts (0.07 < z < 0.1) tend to arrive at a 
normalization that is a factor of ~ 1.5 higher, though 
they state general consistency with lower redshift mea- 
surements given possi ble systematics due to differences 
in methodology, e t c. ([Bell et all l200l lEke et al.l 120051: 
iSmith et all l2009t iHill et al.l 1201 OD . Taken together, 
however, the aforementioned studies appear to point to 
the possibility of an increasing LF normalization from 
z~0.05 toz~0.1. 

Ideally, one would like to study the total luminosity 
density (and hence stellar mass density) as a function 
of distance, rather than just simply the normalization of 
the LF. However, any apparent magnitude limited survey 



at low redshifts is limited on the bright end of the LF 
by poor counting statistics and on the faint end by the 
relatively small(er) volume over which faint galaxies may 
be sampled. These two effects, along with differences in 
methodology, combine to yield the result that the overall 
shape of the NIR LF can be quite different from one 
study to the next. These differences in LF shape can 
lead to substantial differences in integrated luminosity 
density. 

The peak contribution of the LF to the total luminosity 
density occurs at i ~ L*, where all studies feature the 
best statistics. The Schechter (1976) function paramters 
for the normalization (0*) and characteristic luminosity 
(L*) are correlated. Thus, in this study we focus on a 
comparison between our study and those from the litera- 
ture, of the product of the normalization and character- 
istic luminosity {(j)* L*), which amounts to a comparison 
of the peak of the luminosity density distribution across 
studies. 

In this paper, we probe the NIR LF just beyond the 
local volume at redshifts of 0.1 < z < 0.3. To study 
the NIR LF at these relatively low redshifts, high spec- 
troscopic completeness is essential to overcome possible 
biases and to deal with the fact that the errors in the 
photometric redshifts are of the order of the redshifts 
themselves (i.e. ctz/z ~ 1). Our spectroscopic sur- 
vey of i?— band selected galaxies (TJab < 18) is > 90% 
complete over six widely separated fields covering a to- 
tal of 2 deg^ on the sky. Rest-frame band light 
is a good tracer of stellar mass and hence a galaxy's 
-ff— band luminosity (in solar units) is appro ximately 
equal to its stellar mass (in solar units) (Ide Jong 



Tg M IBell fc de Jong|l200TI: IBeh et al.ll2003l : iKhbvetaT 
20081). At low redshifts, Xobserved ~ Arest, and, in the 
NIR, j^T— cor rections are small and nearly independent of 
galaxy type (jMannucci et al.ll2001l ). so this is effectively 
a mass-selected sample. 

Our survey is unique among extragalactic redshift sur- 
veys because of its relative depth, width, and spec- 
troscopic completeness. Deep redshift surveys, such 
as the Cosmic Evolution Survey (C OSMOS / ZCOS- 
MOS, IScoville et all I2007ht ILiUv et al. 2009), the All- 
wavelength Extended Groth strip International Sur- 
vey plus the Deep Evolutionary Exp loratory Survey 2 
(AEGIS + DEEP2. lDa^is et al.ll2007h. and the VIMOS 
VLT Deep Survey (VVDS. ILe Fevre et al.ll2005h all fea- 
ture high-quality NIR photometry and ~ 10, 000 pub- 
licly available redshifts each. However, the spectro- 
scopic completeness at relatively bright NIR magnitudes 
(TJab < 18) is less than 40% in all cases. The Galaxy 
And Mass Assembly survey (GAMA, IDriver et al.llMl 
will feature high completeness to K ^ 18, although the 
current public data release is only complete to K 16. 

The structure of this paper is as follows: We present 
our observations, data reduction, and redshift determi- 
nation methods in Section [51 We discuss our methods 
of spectral energy distribution (SED) fitting to obtain 
photometric redshifts for targets lacking spectroscopy in 
Section |3l We present our results regarding the observed 
NIR LFs in Section H) We analyze the effects of cosmic 
variance on our study and the last decade of studies of the 
NIR LF in the literature in Section [5] We summarize in 
Section |6l Unless otherwise noted, all magnitudes given 
in this paper are in the AB magnitude system (rriAB = 
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23.9 — 2.5 logj^o fi, with in units of /iJy). We assume 
a cosmology of = 0.27, Ha = 0.73, and h = 0.7 in 
our conversion of redshifts to distances. 

2. OBSERVATIONS AND DATA REDUCTION 

Our primary photometry covers six fields in the 
J, H, and if— bands. In Table 1 we show the location 
of our observed fields in right ascension and declination, 
as well as Galactic and supergalactic coordinates. Also 
shown in Table 1 is the area on the sky, number of target 
galaxies, and spectroscopic completeness for each field. 

We made all the H and J— band observations (except 
the CDF-N J— band) with the Ultra Low Background 
Camera (ULBCam) on the UH 2.2 m telescope. We ob- 
served the CDF-N in the JsTs— band with the Widefield In- 
frared Camera (WIRCam) at the Canada France Hawaii 
Telescope 3.6 m (CFHT). CDF-N J-band observations 
with WIRCam were obtained by a group led by Lihwai 
Lin in 2006A. Our group reduced these public data in 
2008. We observed five of our six fields in the if— band 
with WFCam on UKIRT 

We presented the observations, data reduction, star- 
galaxy separation, and bright galaxy c ounts for our 
deep w ide-field NIR imaging campaign in iKeenan et al.l 
(|2010ai) . In that paper we used the bright NIR galaxy 
counts from our survey, in combination with other data 
from the literature, to explore local large-scale structure 
via the slope of the gala xy counts curv e as a fu nction of 
position on the sky. In IKeenan et al.l ()2010bf ) we inte- 
grated galaxy counts from our NIR survey in combina- 
tion with deeper data from the Multi Object InfraRed 
Camera and Spectrograph (MOIRCS) instrument on the 
Subaru Telescope, and other studies from the literature, 
to obtain the best current estimate of the total light from 
galaxies in the NIR. 

The portion of the survey used for this work includes 
only areas where all three bands have uniform coverage 
and consists of approximately 2 deg^ reaching a 5 cr lim- 
iting magnitude of JHK ^ 22 — 23 over ~ 1.8 deg^ with 
another - 0.2 deg^ to JHK - 24. Thus, the NIR pho- 
tometry used for the present work has very high signal 
to noise. 

Before performing photometry on galaxies, we cali- 
brated our NIR images to 2MASS fluxes using rela- 
tively bright point sources in the magnitude range 14 < 
JHK < 16, where 2MASS reports better than 10 a sig- 
nal to noise and our data are well below saturation levels. 
We found our final calibrated images to be quite fiat and 
free from systematic variations i n fiux as a function of po- 
sition. We used the SExtractor ()Bertin fc Arnoutill996D 
MAG_AUTO aperture to do photometry for galaxies in 
this study. We found that this aperture did an excellent 
job of extracting total magnitudes for high signal to noise 
galaxies (such as those used in this study). For a detailed 
description of the data reduction and calibration proce - 
dures used, we refer the reader to IKeenan et al.l ()2010aD . 

Here we selected all galaxies brighter than 18*^ mag- 
nitude in the H—hand in our fields and targeted them 
in a campaign of spectroscopic follow-up. As noted 
in Section [TJ this is essentially a mass-selected sam- 
ple. We note that thermally pulsating asymptotic giant 
branch (TP-AGB) stars may dominate the NIR light in 
youn g (0.2 — 2 Gyr) stellar populations ()Maraston et al.l 
l200l . which may be a source of error in this selec- 



tion for galaxies that co ntain a signifi c ant p opulation 
of young stars. However, iZibetti et al.l (|2012[ ) have re- 
cently found that post-starburst galaxies at z ^ 0.2 
show NIR fluxes, relativ e to op tical, that are consistent 
with IBruzual fc Charlotl (j2003D stellar population mod- 
els, which suggests TP-AGB stars may not have as sig- 
nificant of an impact as previously believed. 

2.1. Fields Observed 

Our first two fields are centered on the Chandra Large 
Area Synoptic X-ra y Survey (CLASXS;Xang et al. 2003; 
ISteffen et al.l 120041 ) an d the Chandra L ockman Area 
North Survey CCLANS: iTVouille et al.l[2008. 2009). Each 
of these fields cover ~ 0.5 deg^ in JHK. These fields 
are located in the Lock man Hole region of low Galac- 
tic HI column density (jLockman et al.l I1986[) . Our 
third field covers a '-^ 0.2 deg^ ar ea centered on the 
Chandra Deep Field N orth f CDF-N. [Brandt et al.|[200TI: 
I Alexander et al.ll2003[ ). The CDF-N contains the Great 
Observatories Origins Deep Survey North (GOODS-N; 
145 arcmin^ HST Ad vanced Camera for Surveys observa- 
tion, [G^^lisco_eLin[2Q03). Our fourth field is the Abell 
370 (A370) cluster and surrounding area (^ 0.4 deg^). 
A370 is a cluster of richness at a redshift of z = 0.37. 
Our fifth and sixth fields are ~ 0.2 deg^ each centered on 
the "Small-Survey- Area 13" (SSA13) and "Small-Survey- 
A rea 17" (SSA17) fr om the Hawaii Deep Fields described 
in ILillv et all (Il99l . 

2.2. Ancillary Photometry 

All the fields included in this work (except A370) are 
covered by the SDSS, which allows us to use uniform 
photometry for fitting SEDs. In add i tion to the fields 
that were presented in IKeenan et al.l (|2010af ) (CLANS, 
CLASXS, CDF-N, A370, and SSA13), we added the 
blank field "SSA17". The NIR data reduction and pho- 
tometry methods f or the SSA17 field we re identical to 
those described in IKeenan et al.l ()2010af ) for the other 
fields, except SSA17 does not include if —band photom- 
etry. 

In the case of the A370 field, we use publicly available 
optical photometry from the Canada France Hawaii Tele- 
scope (CFHT) MegaCam in the u, g, r, and i bands (there 
were no z— band observations for A370). We downloaded 
the catalogs from the MegaCam image stacking pipeline 
website (Megapip43)- The MegaCam u,g,r, and i filters 
are quite similar to those of the SDS^, and small magni- 
tude offsets provided on the Megapipe website bring the 
data into accordance with SDSS magnitudes. 

Catalogs generated through Megapipe are processed 
in the following way: Raw images are reduced or 
"detrended" through the Elixii^ system, which in- 
cludes treatment for dark and bias subtraction, flat- 
fielding, defringing, and basic astrometric and photo- 
metric calibration. Elixir processed images are re- 
trieved and checked for quality before stacking. More 
refined processing for astrometric and photometric cal- 
ibrations are performed and then the images are 
coadded. Catalogs are generated through the appli- 
cation of the source extraction software SExtractor 

® http:/ /cadcwww.dao.nrc.ca/megapipe 

^ http: / / cadcwww.dao.nrc.ca/megapipe / docs / filters.html 

* http:/ /www. cfht.hawaii.edu/Instruments/Elixir 
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TABLE 1 

Coordinates, areas, and number of targets (spectroscopic redshifts / total) 



Field 


CDF-N 


CLANS 


CLASXS 


SSA13 


SSA17 


A370 


R.A.(hh:mm:ss) 


12:36:55 


10:46:54 


10:34:58 


13:12:16 


17:06:31 


02:39:53 


Dec (dd;mm:ss) 


62:14:19 


59:08:26 


57:52:22 


42:41:24 


43:55:44 


-01:34:37 


Galactic 1 (deg) 


125.9 


148.2 


151.5 


109.1 


68.9 


173 


Galactic b (deg) 


54.8 


51.4 


51.0 


73.8 


42.0 


-53.5 


Supergalactic 1 (deg) 


54.7 


52.2 


52.3 


75.3 


65.1 


302.3 


Supergalactic b (deg) 


11.7 


-1.8 


-3.9 


13.6 


49.5 


-25.7 


Area (deg^) 


0.2 


0.5 


0.5 


0.2 


0.2 


0.4 


Spectroscopic Redshifts 


92 


297 


278 


121 


161 


367 


Total Targets {H < 18) 


100 


313 


298 


124 


188 


413 



(jBertin &: Arnouti 119961 ). For a detailed description of 
all processes in the Megapipe reduction, please refer to 
the website (http://cadcwww.dao.nrc.ca/niegapipe). 

Three of our helds (CLANS, CLASXS, and CDF- 
N) are covered in part by Spitzer InfraRed Array 
Camera (IRAC 3.6 — 8 ^m, iFazio et all [200l cam- 
paigns. In the CLANS and CLASXS fields, we use the 
Spitzer Wide-area InfraRed Extragalactic Survey cata- 
logs (SWIRE, Lonsdale et al. 2003) to supplement the 
photometry for 516 of our target galaxies. In the CDF- 
N, we use the IRAC catalogs of iWang etlll (|2010i ) to 
supplement the photometry for an additional 33 targets. 
As a result, almost 40% of our target galaxies have high 
quality photometry in the mid-infrared. 

2.3. Spectroscopic Observations and Reductions 

We observed our target galaxies between 2006 and 
2010 with the Hydra fiber spectrograph instrument 
(jBarden fc Armandrofil[T"995f ) on the Wisconsin-Indiana- 
Yale-NOAO (WIYN) telescope at Kitt Peak in Arizona. 
Hydra is a multi-fiber spectrograph with ^ 90 function- 
ing fibers in the configuration we used. The field of view 
of this instrument is 1 degree in diameter on the sky and 
fibers may be placed a minimum of 37" apart. We con- 
figured the fibers using the whydra configuration code, 
which allows the user to define a minimum number of sky 
positions and guide stars required and allows for targets 
to be weighted by priority. We also used the hydrasim 
simulator code to check that the output of whydra was 
compatible with the software at the telescope. The hydra 
"red" fibers are 2" in diameter, which was ideally suited 
to getting the majority of light from our target galaxies 
onto the spectrograph with a minimum of contamination 
from the sky and neighboring sources. 

For this survey we selected galaxies from our iJ— band 
photometry {H < 18). We masked areas of our images 
near bright objects { Rmncmitudp. < 14.5), as described in 
iKeenan et al.l (|2010a[ ) , to avoid targets that were blended 
with bright neighbors. This selection yielded 1436 target 
galaxies in the range 14.5 < H < 18. 

We configured the spectrograph using the "red" fiber 
bundle and the 316@7.0 grating at first order with the 
GG-420 filter to provide a spectral window of ^ 4000 — 
9500 A with a pixel scale of 2.6 A per pixel. Toward 
the edges of our spectral window the spectra were of- 
ten dominated by noise due to the waning sensitivity 
of the spectrograph on the blue end and the spectrum 
becoming sky-dominated in the red. As a result, the 
useable spectral window was ~ 4500 - 9000 A. We ob- 
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Redshift 
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Fig. 1. — (a) //—band absolute magnitude versus redshift for all 
galaxies in our sample. Black points represent the spectroscopic 
redshift sample and red points represent the photometric redshift 
sample, (b) All redshifts in our six fields. The total for all six 
fields is shown as a black histogram and each field is represented 
individually with colored histograms. 



tained the majority of our target spectra since 2008 
when the bench spectrograph at WIYN was upgraded 
(Bershady et al. 2008). The improved sensitivity of the 
upgraded bench spectrograph allowed us to determine 
redshifts for ~ 90 — 100% of target galaxies in a two 
hour exposure given ideal observing conditions. 

Our exposures were typically 20 minutes each, and we 
would stay on target for 1-3 hours total, depending on the 
observing conditions and the targets observed. For cal- 
ibration we took dark frames, bias (zero exposure time) 
frames, dome fiats, and CuAr comparison lamp expo- 
sures. We employed the iraf task dohydr(^ in the reduc- 
tion of our spectra. This task is specifically designed for 
reduction of data from the Hydra spectrograph and in- 
cludes steps for dark and bias subtraction, fiatfielding, 
dispersion calibration, and sky subtraction. 

^ http:/ /iraf.noao.edu/tutorials/dohydra/dohydra.html 
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2.4. Redshift Determination 

We used the iraf ta sk xcsao within the rvsao package 
(jKurtz fc Minklll998f ) to determine redshifts for our ob- 
served spectra. The xcsao task is an interactive tool that 
allows the user to determine redshifts by cross-correlating 
a spectrui n against template spe ctra using methods de- 
scribed in iTonrv fc Davi^ (|1979( ). This task allows the 
user to manually mask portions of the spectrum and cy- 
cle through a variety of candidate fits to find the best 
match. We manually determined a redshift for all of the 
galaxies in our spectroscopic sample using xcsao in this 
manner. 

At the relatively low redshifts of our sample galax- 
ies, the commonly identifiable absorption features in 
our spectra were Call H (3969 A) and K (3934 A), 
the 4000 A break, the Fe G-band feature (- 4100 A), 
H7 (4342 A), H/3 (4865 A), and the Na D lines 
(5890/5896 A). If emission lines were present, our spec- 
tra typically captured the Ha (6563 A) and [Nil] 
(6548/6584 A) complex along with the [SII] dou- 
blet (6719/6730 A) and/or the [Oil] (3727 A), [OIII] 
(4961/5007 A) and (4865 A) emission hues. 

For ~ 50 galaxies we targeted objects that already had 
a cataloged redshift in the SDSS archive. We find 100% 
agreement with the SDSS redshifts to within an rms er- 
ror of 0.0005. We supplement our catalogs with an ad- 
ditional 32 redshifts from the SDSS for other targets in 
our fields. We also supplement our catalogs with 129 
redshifts for objects in our fields that were either previ- 
ously published in the literature or obtained from other 
spectroscopic campaigns within our research group. 

For the CDF-N fiel d, many of th ese supplementary red- 
shifts are cataloged in lBarger et a l. (2008) and references 
therein. In the CLANS and CLASXS fields, man y of the 
supplementa ry redshifts are cataloged in Trouill e" et al.l 
(I2OOI [2OO9I). Manv re dshifts in the SSA13 field are cat- 
aloged inl^wieiLal] ([2004 ) and in the A370 field in I. 
Wold et al. (2012, in preparation). The remainder come 
from an ongoing campaign of followup spectroscopy in 
our fields using the DEep Imaging Multi-Object Spec- 
trograph (DEIMOS) on the Keck telescope (L. Cowie, 
private communication). 

In Figure [Iji, we show band absolute magnitude 
versus redshift, where black points represent the spectro- 
scopic sample and red points represent photometric red- 
shifts (see Section[3]). In Figure[TlD, we show the spectro- 
scopic redshifts in our fields. The black histogram shows 
the total for all six fields, and the colored histograms 
show each field individually. 

For our redshifts obtained with the Hydra instrument 
(1155 total), we assigned a simple confidence flag to the 
determined redshift to indicate whether it was a high- 
confidence redshift (1) or questionable (0). In reality 
there may be some ambiguity, in relatively low signal to 
noise spectra, as to when one can robustly determine a 
redshift versus when a redshift may be determined but 
could be questionable. To keep things simple in this case, 
we assign high confidence to redshifts where a minimum 
of two absorption or emission features could be clearly 
identified in different parts of the spectrum simultane- 
ously. We assign questionable confidence when one or 
more features appear identifiable in a spectrum, such 



that a redshift could be defined, but that other possi- 
ble fits to the spectrum cannot be completely ruled out. 
With this system we assign high-confidence to ~ 92% of 
the redshifts in our Hydra sample. 

3. SPECTRAL ENERGY DISTRIBUTION FITTING USING 
GAZELLE 

We used the GAZELLE code (R. Kotulla et al. 2012, 
in preparation) to fit the SEDs of our target galax- 
ies in order to obtain photometric redshifts for galaxies 
that lack spectroscopy. The GAZELLE code uses the 
GALaxy EVolution (GALEV. EotuUa et al.l 12009) stel- 
lar population synthesis models as the basis for fitting 
SEDs. GALEV is an evolutionary synthesis code that 
allows the user to generate models to describe the evo- 
lution of stellar populations in general, whether they be 
resolved populations (as in the case of star clusters) or 
unresolved populations in galaxies. 

Given a set of initial conditions, GALEV is able to 
track the integrated light properties of galaxies over 
cosmological timescales. Users are now able to gener- 
ate their own models directly via the GALEV website 
(http://www.galev.org). To generate a model, the user 
first chooses a galaxy type (E, Sa, Sb, Sc, Sd) or specifies 
a particular star formation history (e.g., constant, expo- 
nentially declining etc.), including an option to upload a 
user defined star formation history. 

GALEV is unique among evolutionary synthesis codes 
in that it allows the user to select a "chemically consis- 
tent" treatment by simultaneously tracking the spectral 
evolution of stellar populations and the chemical evolu- 
tion (metallicity) of the gas component of galaxies. Al- 
ternatively, one can choose from five fixed metallicites in 
the range —1.7 < [Fe/H] < -t-0.3. The user may also vary 
the initial mass function (IMF) , switch on or off the emis- 
sion from gas, choose a Galactic extinction prescription, 
and vary parameters regarding the details of star forma- 
tion efficiency and burst duration, age, e-folding time, 
etc., for starburst models. Lastly, the user chooses a cos- 
mological parameter set, desired outputs of the model, 
and the photometric filter profiles to be applied to model 
SEDs. 

In an effort to fit model SEDs to photometric data, 
a wide selection of models is required to cover a range 
of star formation histories, ages, metallicities, galaxy 
types, and so on. In this work, we chose 83 indepen- 
dent GALEV models to use in fitting our photometry. 
20 of these models are the types available through the 
GALEV website for five galaxy types (E, Sa, Sb, Sc, 
and Sd). Each galaxy type is modeled at three fixed 
metallicities ([Fe/H] = —0.7,-0.3,0) and once in the 
"chemically consistent" evolving metallicity format. The 
other 63 models are starburst models that attempt to 
capture a wide range in post-burst age. The starburst 
models are not publicly available through the GALEV 
website. They are described in detail in the disserta- 
tion of Ralf Kotulla and were provided to us directly 
by him. Altogether, we use burst models covering 21 age 
steps (post-burst 0.002-9 Gyr) at three fixed metallicities 
([Fe/H] = -0.7,-0.3,0). 

GALEV does not yet handle emission from dust, so 
that in low-redshift samples such as this work, the longer 
wavelength Spitzer bands contribute relatively little to 
the SED fits. This is because the GAZELLE code ignores 
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photometry in bands where dust emission is expected to 
dominate the signal (> 3 ^m). 

GAZELLE uses a algorithm to fit models to the 
observed data. The values for each model are then 
transformed into normalized probabilities. The model 
with the highest probability is the one chosen for com- 
puting redshift, mass, metallicity, and so on. The formal 
1 a uncertainties in redshift, mass, metallicity etc. are 
derived by finding the minimum and maximum values 
associated with the models in the top 68% of the nor- 
malized probability distribution. 

GAZELLE includes an option to run with a user de- 
fined redshift for galaxies. This mode allows the user 
to "calibrate" GAZELLE given the difference between 
observed photometry and the output model photometry 
from the best-fit SED. In this mode, we compared the 
output model magnitudes from GAZELLE with the in- 
put observed magnitudes and found very good agreement 
in the optical {\mobserved - rrimodeil < 0.03), but system- 
atic offsets of ~ 0.1 — 0.2 in the NIR. We then included 
these offsets in the GAZELLE input parameter file to 
account for these systematics, such that, when we re-ran 
the code, the model output magnitudes agreed with the 
input magnitudes in all bands with systematics of < 0.03 
magnitudes. These adjustments served to improve the 
spectroscopic versus photometric redshift relation. 

At the relatively low redshifts of our sample, the er- 
ror in photometric redshifts can be of the same order 
as the redshifts themselves. Despite these relatively 
large uncertainties, we find that the photometric red- 
shifts derived for galaxies that lack spectroscopy fall 
within the range of expected values from the spectro- 
scopic sample. We note that > 50% of our targets 
have a \zspec — Zphot\/{^ -I- z) < 0.05 and > 80% have 
\zspec — Zphot\/{^ + z) < 0.1, which could be considered 
fairly typical for determining photometric redshifts. 

4. THE NIR LUMINOSITY FUNCTION 

The NIR LF of galaxies has been relatively well estab- 
lished at low redshifts of (z) ^ 0.025- 0.06 (ICole et al.l 
[2OOII: IKochanek et~aLl [2OOII: IJones et al.l I2006D . At 
slightly higher redshifts of (z) ~ 0.06 — 0.1, large sur- 
veys have also provided constraints on the LF and tend 
to arrive at a normalization that is a factor of ~ 1-5 
times higher than lower redshift studies (iBell et al.|[2003l : 
lEke et al.ll2005l: ISmith et"alll200l IhTU et al.ll2010D . 

At the mean redshift of our sample (z ~ 0.2) the LF is 
not particularly well constrained because, in this redshift 
range, depth and wide area are simultaneously required 
to sample the entire LF. Furthermore, high spectroscopic 
completeness is essential at these redshifts to overcome 
possible biases, and to avoid the relative uncertainty in 
photometric redshifts {<Tz/z ^ 1), which lead to large 
uncertainties in absolute luminosities. 

Here we combine our fields to study the NIR LF nor- 
malization at 0.1 < z < 0.3 and to compare with lower 
redshift values. As noted above, our sample is selected 
in the iJ— band, but we have photometry for the en- 
ti re sample in the J a nd iiT— bands as well (presented 
in lKeenan et al.|[2010af) . Thus, we are able to derive LFs 
in J, H and K to compare with a wide range of selections 
from the literature. 

First, we derived LFs for the combination of all six of 
our fields in the iJ— band using four different estimators 



for both the LF and its normalization. From these meth- 
ods we devise a hybrid scheme of deriving the LF that 
combines two different LF estimators. We compute the 
J, H, and if —band LFs in this way to compare the over- 
all shape and normalization with other studies from the 
literature. In the range 0.1 < z < 0.3, we have 812 galax- 
ies. This subsample is 92% spectroscopically complete. 
The only other study to focus specifically o n the NIR LF 
in thi s redshift range was performed by iFeulner et al.l 
(|2003f ). Our sample represents a factor of four increase 
in the number of galaxies used in their work. 

4.1. Methods for Calculating the Luminosity Function 

Several methods exist for calculating the LF of galax- 
ies. The most co mmonly as s umed functional form of the 
LF is that of the iSchechteil (|1976[ ) function, which takes 
the following form 



$(L)dL 



exp 



L\ dL 



(1) 



which may be converted to be in terms of absolute mag- 
nitudes 



±l_ ^ -[^Q-0.4(Af-M*) 

L^ 



(2) 



such that 



$(M) = 0.41n(10)(/)* 



-[^q0.4(M*-M) 



(a+l) 



exp(100-4(M--M)) 



(3) 



The Schechter function parameter L^, (or M*) repre- 
sents the luminosity of galaxies at the knee of the LF, 
while (j)* determines the number density of galaxies, 
and a is the faint-end slope. Methods for deriving the 
LF include parametric and non-parametric forms. In this 
work we use a hybrid combination of parametric and non- 
parametric estimators to determine the Schechter func- 
tion parameters, but first we compare four methods of 
LF estimation independently. 

4.2. Parametric LF Estimation 

The most commonly used parametric method is that 
of iSandage et al.l (|1979D (STY). This method employs a 
maximum likelihood estimator (MLE) and the assump- 
tion that the LF takes the Schechter function form. An 
advantage of this method is that no binning is required. 
Instead, a probability density, p{Mi, Zi), is calculated for 
each galaxy with absolute magnitude Mi at redshift Zj, 
where Mfaint(2i) is the faintest galaxy visible in an ap- 
parent magnitude limited survey at redshift Zi , such that 



$(M,) 



X 



^{M')dM' 



(4) 



The likelihood function is then the product of the prob- 
ability densities for all galaxies in the survey 



C = \{p{M,,z, 



(5) 
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This product is maximized to obtain the best fit 
Schechter function parameters. In this method, the den- 
sity functions cancel in the calculation of the probability 
densities. Thus, this technique is insensitive to density 
inhomogeneities, and the normalization (0*) must be de- 
termined via some other method. 

4.3. Non-parametric LF Estimation 

Non-parametric LF estimators include binned and non- 
binned methods, and no a priori assumption of the 
LF functional form is required. However, for purposes 
of comparison with other studies and other methods a 
Schechter function is usually fit to the resulting LF esti- 
mation. 

The sim plest non - param etric approach is the 1/Vmax 
method of iSchmidO (|1968f l. This method assumes a ho- 
mogeneous distribution of galaxies, such that each galaxy 
is assigned a weight (l/Vmax) based on the fraction 
(Vmax) of the total survey volume in which that galaxy 
can be observed, given its absolute magnitude. The sum 
of the weights for galaxies in an absolute magnitude bin 
is then the number density, $(M), for that bin. This 
method has the advantage that the normalization of the 
LF is determined directly, but it has been shown that 
this method may underestimate the number density of 
galax ies near the flux limit of a sample (Page & Carrera! 

Another non - param etric approach is the method 
of lLvnden-Belil (flOTl . The method does not require 
binning, and it does not require the assumption of a ho- 
mogeneous distribution of objects. Thus, as for STY, 
the normalization must be determined by other means. 
In the method, the cumulative luminosity function 
is given by 



/ ^{M)dM = n 

J —oo 



(6) 



In a sample with a redshift lower limit of ^min, the 
value for the parameter C^T for each galaxy at (Mfe, Zfe) is 
the total number of galaxies that are both brighter than 
Mk in absolute magnitude and within the redshift range 
(^min < Zk < Zk,Hmit), whcrC Zk, limit is the maximum 
redshift that a galaxy with absolute magnitude may 
be observed in the survey. In other words, 

i 

where the window function W{Mk) is defined as 



W{Mk) 



1, M <Mk and z„ 
0, otherwise 



< Zk < Zk, limit 



The resulting cumulative luminosity function may then 
be binned and fit with a Schechter function for compar- 
ison with other studies. 

Another popular non-parametric approach is the 
"Step- Wise Maximum Likelihood" estimator (SWML, 
lEfstathiou et al.l 119881) . This method requires binning 
but no a priori assumption of homogeneity. In this case, 



0) 
Q 

0) 
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-24 -23 -22 -21 -20 -19 

Mh - 5log,o(h) 

Fig. 2. — A comparison of the band LFs derived for our sam- 
ple using four different methods. The only method that establishes 
the normalization is the 1/Vmax method, so all four methods have 
been set to that normalization for comparison. 



the LF is parameterized into Np steps of width AM such 
that 



$(M) = $ft 



l,...,iV„ 



(9) 



Two window functions, W{x) and H{x), are defined 
such that 



Wix) 



1, - AM/2 <x < AM/2 
0, otherwise 



(10) 



Hix) 



1, x> AM/2 

{x/AM + 1/2), -AM/2 <x < AM/2 
0, x<-AM/2 



(11) 

Maximizing the l ikelihood function (please see 
lEfstathiou et al.lll988l for details) then yields the follow- 
ing form of the differential LF 



j:^W{Mk-M,) 



-H"[Mii„(z.)-Mfc] 



(12) 



This estimator is then applied iteratively, each time re- 
placing the input values for $j with the output $fe values 
from the previous iteration until sufficient convergence is 
achieved. 

4.4. Comparison of Four LF Estimators 

In Figures [2] and [HI we compare the LFs and Schechter 
function parameters AI* and a determined from the four 
LF estimators described in Sections |4. 21 and r4. 31 (1/1/^^^. 
STY, SWML, and C") . 

Figure [5] shows a comparison of the iJ— band LF de- 
rived for our sample using the four different methods. 
The l/Vmax method is shown in black, the STY method 
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Fig. 3. — A comparison of the iJ— band LF parameters M* and 
a derived for our sample using four different methods. The elHpses 
show the 1 cr errors in these parameters for each method. 



in red, the C~ method in blue, and the SWML method 
in orange. The only method that establishes the normal- 
ization is the 1/Vmax method, so all four methods have 
been set to that normalization for comparison. 

Figure [3] shows a comparison of the iJ— band LF pa- 
rameters M* and a derived for our sample using the four 
different methods. The same color designations apply as 
for Figure El 

We expect the a parameter for our work to be close to 
the value of a « — 1 found in large low-redshift studies, 
as evolution in the LF should be minimal since z ~ 0.2. 
Given our comparison of different methods above, we 
decided to use the STY method to determine the a pa- 
rameter in the sections that follow. We determine the 
error in our calculation of a by fitting a 2-dimensional 
Gaussian to the distribution of values for the likelihood 
function in a, M* space. We find that a 2-d Gaussian 
is a good fit to the distribution. We show the 1 tr error 
ellipse of this Gaussian is that shown in Figure [3] 

In further exploration of fitting a with STY, we found 
that the method is sensitive to the input apparent mag- 
nitude limit of the survey. Our sample is i?— band se- 
lected, such that the magnitude limit is well known in 
the i?— band but less well defined in the J and A'— bands. 
Thus, we first determine a in the i/— band (a = 0.91 ± 
0.11) and then hold that value fixed in our determination 
of M* using STY for J, H and K. 

In Section [SJ we discuss a comparison with the COS- 
MOS cone mock catalogs to study cosmic variance. We 
use the mock catalogs to test our method of determining 
a and M* with STY. To do this we combined all the 
mock catalogs to form a simulated survey 24 times the 
size of our survey. We applied the above methods to de- 
termine a and M*. We then created mock surveys of 6 
random sightlines the same size as our survey and per- 
formed the same process to test whether the values for a 
and M* recovered from the smaller survey were consis- 
tent with the global values. We found that the recovered 



values for a and M* for the smaller mock surveys were 
in 1 (T agreement with the global values determined from 
combining all the mock catalogs. 

4.5. Determining The LF Normalization 

Of the four methods discussed above for determin- 
ing the LF Schechter function parameters, the l/Vmax 
method is the only one that gives the normalization di- 
rectly. The other methods require that the normalization 
be determined by other means. 

Several methods have been developed for determining 
th e LF normaliza tion, and these are studied in detail 
by iWillmed ()1997[ ) via simulated data. They find that 
the various density estimators yield roughly equivalent 
results, and all tend to underestimate the true density of 
objects by some ^ 20%. The parameter (j)* is related to 
the mean density (n) of the sample through 



(13) 



where M^right and Mfaint are the absolute magnitude 
limits of the survey. 

One commonly used method for determining the 
mean density is the minimum variance estimator of 
IDavis fc Huctral ([l982h 



^dVS{z)w{z) 



(14) 



where Ngai is the number of galaxies in the sample, S{z) 
is the selection function, and w{z) is a weighting func- 
tion for each galaxy defined by the inverse of the second 
moment (J3) of the two-point correlation function, ^(r). 
The selection function, S{z)^ is defined as 



S{z) = 



(15) 



where M^axizi) is the faintest absolute magnitude de- 
tectable at redshift z^. 
The weighting function is defined as 



w{z) = 



1 



l + nJaSiz) 



Jo 



(16) 



In the case that all weights are set to a value of 
1, the minimum variance estimator described above 
reduces to the so called "^3" estimator proposed by 
IDavis fc Huchral (|I982. ) 



Nt 



(17) 



Ano ther estimator proposed by IDavis fc Huchral 
Ip8l is 



ni 



S{z) 



dz 







]/ 






Jo dz 



(18) 



which in the c ase of 1 galaxy pe r bin reduces to the form 
used bv ,Efstathiou et all (| 19881) : 
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TABLE 2 

BAND NIR LF NORMALIZATION </)* 

DERIVED USING FOUR DIFFERENT 
i\/r 17 'runny t) 



Method 


(f>* X 10=^ 


1/ Vmax 


15.3 ± 1.6 


^minvar 


15.6 ± 1.4 


ni 


14.9 ± 1.5 




15.7 ± 1.3 



cthods listed 



plu 

EE 



three methods (ti^ 



of Mpc 
are 1/V™ 
L. "3) 



(19) 



We show a comparison of the normalization derived 
using the aforemention ed mean density estimators of 
iDavis fc Huchral (jl982) alongside the result for the 
^/Vmax method in Tabic [2l In all cases, we have de- 
termined a and M* using the STY method. We find 
that all four estimates of the LF normalization agree to 
within 1 a statistical errors. 

4.6. A Hybrid Method for Determining the NIR LF 

All LFs for our observed data presented in the following 
sections are calculated using the method described in 
Section [4.4[ where, using STY, we first determine a from 
the iJ— band sample and then M* for each bandpass. We 
then determine 4>* by fitting a Schechter function to the 
1/Knax results in each bandpass with a and M* fixed at 
the values determined with the STY method. 

The ^/Vmax method requires binning. Thus, we com- 
pared results from bin sizes ranging from 0.2 — 1 magni- 
tudes. We found the results to be consistent to within 1 a 
(statistical uncertainty in the resulting Schechter func- 
tion fit). We settled on a bin size of 0.4 magnitudes to 
provide the best counting statistics over a wide range in 
absolute magnitudes. 

The STY method does not allow for any test of good- 
ness of fit between data and model. We fit the Schechter 
function parameter (jf to the 1/Vmax binned data, with 
a and M* fix ed, using th e mpfit .pro least-squares fitting 
procedure of i Markwardtl (|2009| ). which is based on the 
MINPACK^^ algorithm of Jorge More. The resulting 
fits yield a. ^ 15.5 for 10 degrees of freedom, suggest- 
ing that the Schechter function is a good fit to the data 
at the ~ 90% confidence level. In Section [5l we discuss 
the systematic uncertainties in our determination of the 
LF normalization due to cosmic variance. 

To calculate absolute magnitudes we modified the 
observed apparent magnitude by a distance modulus 
{DM), a iC— correction K{z), and an evolution correc- 
tion E{z), such that M = m - DM{z) - K{z) + E{z). 

4.7. NIR K — corrections K{z) 

At low redshifts, i^T— corrections are small and 
nearly independent of galaxy type in the NIR 



• J (This Study) 
0.018 - ,H (This Study) 
- (This Study) 



iff 



Cole+ 2001 J O 
Cdle+ 2001 K O 
Bell+ 2003 K □ 
Johes+ 2006 J o 
Jdnes+ 2006 H 
Johes+ 2006 K o 
Smith+ 2009 K x 
Hill+ 2010 J A 
Hill+ 2010 H 
Hill+ 2010 K A 
Hudng+ 2003 K S> 
Eice+ 2005 J -Ci ' 
Eke+ 2005 K 
Kdchdnek+ 2001 K + 
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-21.0 -21.2 -21.4 -21.6 -21.8 -22.0 -22.2 
M* - 5Log,„(h) 

Fig. 4. — 0* versus M* for various studies. Studies in the 
J— band are shown in blue, /f— band in green, and if— band in 
red. The error bars show the 1 a statistical errors in each study. 
We note the trend toward a lower value for 0* at brighter values 
for M* . 



CMannucci et"aI1[200l. U sing SDSS and UKIDSS data, 
Chilingari an et aP ()2010( ) showed that, at low redshifts 
(z < 0.5), accurate ii'— corrections may be calculated us- 
ing low order polynomials and inputs of only redshift and 
one observed NIR color. They provide a ii'— correction 
calculatoiF*]. which we used to compute if —corrections 
for galaxies in our sample given redshifts and NIR colors. 

We find the calculator output values to be in agreement 
with the expected NIR j^T— corrections for gal axies at 
these redshifts derived by iMannucci et ap (|2001[ ) . Given 
the analysis of iChilingarian et al.l ( 2O1O0 in comparing 
their calculator outputs with more rigorous SED fit- 
ting methods for determining X— corrections, we expect 
our magnitude errors associated with the if- corrections 
should be < ±0.1. 

4.8. Evolution Corrections E{z) 

Evolution in the rest-frame NIR for galaxies at rela- 
tively low redshifts is expected to be weak but still has 
a significant effect on the normalization of the LFs, par- 
ticularly at z > 0.1. The sim plest form of the ev olution 
correction is E{z) = Qz f e.g..lSmith et alll2009fl . where 
Q is a positive constant. iBlanton et al.l (|2003 t) showed 
that in the NIR, Q = 1 agrees well with stellar popu- 
lation synthesis models. Thus, for this work we adopt 
Qj = Qh — Qk — 1, such that E{z) = z in all bands. 



4.9. 



The Correlation Between Schechter Function 
Parameters 



The Schechter function parameters (p* and M* are cor- 
related, as shown in Figure ID Thus, rather than a direct 
comparison of (f>* values from various studies, we consider 
the product of the normalization and the characteristic 
luminosity {(j)*L*) in the following discussions of compar- 
isons between our study and those from the literature. 

The Schechter function parameter a is also correlated 
with (p* and M* but uncorrelated with the product of the 



http:/ /www. netlib.org/minpack 
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Fig. 5. — The normalized product (p*L* versus a for various 
studies. Studies in the J— band are shown in blue, band in 
green, and X— band in red. The error bars show the 1 a statistical 
errors in each study 



normalization and characteristic luminosity, as shown in 
Figure [5l Thus, our comparison of the product (jj* L* 
between studies is unaffected by the variations in a from 
one study to another. 

4.10. Comparing Studies With Different Photometric 

Methods 

Here we briefly discuss the photometric methods em- 
ployed in different studies that we compare with in this 
work. For all our photometry, we cali brate to 2MASS 
point sources and use the SExtractor (.Bertin fc ArnoutsI 
I1996D MAG _AUTO aperture, which fit s a Kron-like 
Kronl [T980[ ) ellipse to each galaxy. In iKeenan et al.l 



2010af) we found that this aperture does an excellent 



job of retrieving total magnitudes for galaxies at high 
signal to noise (galaxies in the current st udy are 4 — 5 
magni tudes above our detection limits of iKeenan et al.l 
\2mM1 . 

In the case of Kochan ek et al.l (I2001D. iCole et al.l 
(l200l. IB ell et al. (.2003,) . lEke et al.l (I2005D . and 
iJones etal . (2006), 2MASS Kron or similar apertures 
were used for photometry, and we expect this to be con- 
si stent with our w ork. 

iHill et all (|2010( ) apply Kron aperture p hotometry to 
the U KIDSS Large Area Survey (LAS). iSmith et al.l 
(|2009f) also work with the UKIDSS LAS, but they use 
Petrosian magnitudes. They compare the difference be- 
tween Petrosian and Kron magnitudes for ^ 7, 000 galax- 
ies that also have photometry in 2MASS. They find that 
2MASS Kron magnitudes are, in general, ~ 0.2 mag- 
nitudes brighter than their UKIDSS Petros ian counter- 
arts. We adjust the value for M* from ISmith et al.l 
20091) to accommoda te this difference. 
'Hu ang et al.l (|2003l ) use circular 8" apertures corrected 
to 20" apertures for galaxies at i^vcga > 13. While it 
is difficult to assess exactly how this method compares 
to our photometry, the errors on their measurement are 
large and so this d oes not affect our an alysis. Similarly, 
the sample size of iFeulner et al.l (|2003[ ) is small enough 
that the associated uncertainties span essentially the en- 



tire range of measurements of (f>* L* from the literature. 

4.11. The Effects of Magnitude Errors on the LFs 

The fundamental selection for this work comes from 
H—hand photometry we obtained with the Ultra-Low 
Background Camera (ULBCam, Hall et al. 2004) on the 
University of Hawaii 2.2 m telescope. Our final 7?— band 
mosaics are highly uniform and 4 — 5 magnitudes deeper 
than the selection limit for this work {H < 18). We 
calibrated our photometry for all six fields using high 
signal to noise point sources in 2MASS and found no 
discernable systematics and a random magnitude er- 
ror of < 0.05. For the galaxies in this survey, which 
are extended and somewhat fainter than the calibration 
sources, our random magnitude error is ~ 0.1. 

4.12. The Effects of Redshift Errors on the LFs 

Our redshift sample consists of ~ 8% photometric red- 
shifts and ~ 8% questionable spectroscopic redshifts. As 
demonstrated by comparison with the SDSS, the high- 
confidence spectroscopic sample could be considered ac- 
curate to within Az ^ 0.0005. Thus, here we consider 
only the error associated with questionable and photo- 
metric redshifts. 

If we assume, as an upper limit, that all questionable 
and photometric redshifts have an associated random er- 
ror of Az/ (1 -|- z) ~ 0.1, then this would imply a random 
error in a calculated absolute magnitude of ^ 1—1.5. 
Although a large effect for any individual galaxy, the 
contribution of this effect to the mean magnitude error 
in any given magnitude bin would be of the same order 
as the contribution from the random magnitude errors 
associated with the spectroscopic sample. 

4.13. The Effects of Evolution Corrections on the LFs 

ICole et al.l (|2001D showed that for NIR LFs at (z) ~ 
0.05, the evolution correction has the effect of raising the 
normalization by ^ 10% and dimming M* by ~ 0.1 mag- 
nitudes (thou gh these ch anges were of the same order as 
their errors). IJones et~a l. (2006) did not apply an evolu- 
tion correction for their sample at (z) ~ 0.05, and they 
measure a significantly low er normalization and brighter 
M* than ICole etall (|200l . In general, however, apply- 
ing an evolution correction does not significantly change 
the product of the normalization and characteristic lu- 
minosity ((j)*L*). 

4.14. Effects of Magnitude Errors in General 

The E{z) and K{z) errors, along with the aforemen- 
tioned magnitude errors, amount to a total error in the 
calculated absolute magnitudes of ~ 0.3. If these errors 
are truly random, then the effect on the LFs will be to 
shift galaxies from one bin to another, but the random- 
ness of the process should essentially cancel itself out. 

4.15. The NIR Luminosity Function at 0.1 < z < 0.3 

In Figure [6l we show the J, H, and — band LFs of 
galaxies in our six fields combined for the redshift range 
0.1 < z < 0.3. In the left-hand panels (a, c, and e), the 
entire LFs are shown to demonstrate the overall shape 
compared with lower redshift selections from the litera- 
ture. In the right-hand panels, we zoom in on the nor- 
malization (M* ~ —21.7) to highlight the differences be- 
tween this study and other studies from the literature. 
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Fig. 6. — The J,H, and X— band LFs of galaxies in our six fields combined for the redshift range 0.1 < z < 0.3. We derived the 
parameters M* and a via the STY method and then applied the 1/Vmax method and fit a Schechter function with M* and a fixed to 
derive the normalization 0*. In the left hand panels (a, c, and e) we plot the entire LFs to show overall shape for comparison with lower 
redshift LFs from the literature. The colored lines show LFs from the studies denoted in the plots. In the right hand panels, we zoom in 
on the normalization to highlight the differences between the normalization in this study and selections from the literature. The solid lines 
show the Schechter fits to the combined dat a for all s ix fields (five fields in the /S'— band where SSA17 lacks photometry). The error bars 
shown are the 1 a Poisson errors taken from [Gehrelsl (IIQSS) . The normalizations and M* values to the fitted Schechter functions for the 
combination of all six fields are given in Table |3l 
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Fig. 7. — The product of the Schechter function normalization and characteristic luminosity {<f>* L*) for NIR LFs in this study (filled 
circles) and selections from the literature. We have employed a color correction, such that all data from J, H, and K ap pear vertically on 
the same scale. We then normalized the if>* L* values to an error weighted average of the large studies at (2) = 0.057 of , Cole et al.l 120011 ) 
and lJones et al.l lj2006). All measurements from our study and from the literature for the J— band are shown in blue, band in green, 
and X— band in red. In each case, (p* L* is plotted at the median redshift for galaxies in a given study. Where studies (including this one) 
have covered more than one bandpass, the J— band <j)*L* is plotted at the median redshift of the study and H and/or ii"— band values are 
shifted slightly along the abscissa for clarity. The error bars show the 1 cr statistical errors in the product <f>* L* . The values for ip* and 
L*(Af*) used to generate this figure are given in Table|3] The solid line shows the void radial density profile from lBo leiko fc Sussman (20TJ) 
that they claim would be sufficient to provide an apparent acceleration of the expansion of the universe observed via type la supcrnovae. 
IBoleiko SussmanI l|2011l ) presented this void p r ofile a s a density contr ast; (p(r))/po. We have converted from their density contrast to 
(/i*_L* by normalizing to the results of lCole et al.l 120011) and lJones et al.l lf2006. ) . 
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TABLE 3 

NIR LUMINOSITY FUNCTIONS PARAMETERS 0*^ AND Af*'' FOR THIS STUDY AND SELECTIONS FROM THE LITERATURE 



Study 




J- 


-band 


H- 


-band 


K- 


-band 








0* X 10^ 


Al* 


4>* X 10^ 




<t>* X 10^ 


M* 


CV<= 


This Study<i 


0.2 


15.9±1.7 


-21.49±0.05 


15.3±1.6 


-21.67±0.05 


14.6±1.7 


-21.56±0.06 


14% 


Feulner et al. 2003 


0.2 


14.9±2.2 


-21.56±0.24 






11.1±1.2 


-21.95±0.24 


27% 


Huang et al. 2003 


0.14 










13.0±2.0 


-21.86±0.08 


11% 


Smith'' et al. 2009 


0.1 










16.6±0.8 


-21.49±0.05 


2% 


Bell et al. 2003 


0.08 










14.3±0.7 


-21.45±0.05 


10% 


Hill et al. 2010 


0.08 


15.5±2.1 


-21.46±0.13 


14.9±1.5 


-21.89±0.08 


15.6±1.6 


-21.66±0.10 


30% 


Eke et al. 2005 


0.07 


13.9±0.6 


-21.50±0.05 






14.3±0.8 


-21.59±0.04 


6% 


Cole et al. 2001 


0.06 


10.4±1.6 


-21.47±0.02 






10.8±1.6 


-21.60±0.02 


8% 


Jones et al. 2006 


0.05 


7.1±0.1 


-21.96±0.04 


7.2±0.1 


-22.16±0.04 


7.5±0.1 


-21.99±0.03 


5% 


Kochanek et al. 2001 


0.023 










11.6±1.0 


-21.55±0.05 


11% 



^ all 0* values arc given in units of Mpc^"^. For this study wc use a fixed a = -0.91 derived in the iJ-band using the STY method as 

described in Section fm 

^ all A-/* values arc given as M - 51ogiQ(h) 

^ Cosmic variance estimates in percent using the empirical formula of lUriver &; Hobothaml |2U1U 1 

The errors in represent the errors in the Sehcehtcr fits plus the systematics due to cosmic variance. 

The A-/* value listed for Smith et al. is adjusted to be 0.2 magnitudes brighter than their published value to accommodate the difference 
between Petrosian and Kron apertures as described in Section r4. 101 



The solid lines show the Schechter fits to the combined 
data for all six fields (five fields in the X— band where 
SSA17 lacks photometry). 

Given the i?— band value for M* ^ —21.7, we can see 
from Figure [1] that we are sampling M ~ M* galaxies 
out to a redshift of 2; ^ 0.3. We also note from the 
distribution shown in Figure [T] that we don't have a sig- 
nificant number of M ~ M* galaxies below a redshift of 
z ~ 0.1, such that the normalization of the NIR LF is 
best sampled in this survey for 0.1 < z < 0.3. 

In Figure [71 we show the product of the LF normaliza- 
tion and characteristic luminosity {(p* L*) as a function 
of redshift from this study (filled circles) and from some 
literature studies. All measurements in the J— band are 
shown in blue, 7?— band in green, and band in red. 
We applied a color correction so that all three bandpasses 
would appear on the same scale vertically. We then nor- 
malized the (j)* L* values to an err or- weighted average of 
the large studies a t (z) = 0.057 of ICole et al.l ((20011 and 
iJones et"ari ()2006f ). In each case, is plotted at the 
median redshift for galaxies in a given study. Where 
studies (including this one) have covered more than one 
bandpass, the J— band (t)*L* is plotted at the median 
redshift of the study, and H and/or if— band values are 
shifted slightly along the abscissa for clarity. The error 
bars show the 1 a statistical errors in </>* for each study. 

The solid black lin e in Figure [7l shows t he vo id radial 
density profile from 'B oleiko &: SussmanI (|2011l ). which 
they claim could produce the apparent acceleration of 
the expans ion of the universe observ ed via type la su- 
pernovae. iBoleiko fc SussmanI (|2011f) present this void 
profile as a density contrast profile, (p(r))/po- We con- 
vert their density contrast profile to normalized <I)*L* by 
norm alizing to the aver age of the large low-red shift stud- 
ies of ICole et all (|200l and IJones et all (|2Q06t ). We note 
that such a radial density profile is not ruled out by cur- 
rent measurements of the NIR LF normalization. 

We find that our measured value of the product (f)*L* 
is ^ 30% higher than the mean value of studies at 
(z) ~ 0.05. Our measurement could be considered a 
conservative underestimate of the true value of (jf L* for 
0.1 < z < 0.3, because our survey avoids known galaxy 



clusters in this redshift range. 

5. COSMIC VARIANCE 

The relative scatter in measurements of the NIR LF 
normalization over the past ten years could be arising 
from a variety of sources, including differences in fit- 
ting or normalization methods, differences in K{z) or 
E{z) corrections, photometry errors or incomplete spec- 
troscopy, and cosmic variance. Most likely, these differ- 
ences arise from some combination of the aforementioned 
sources, but in the majority of modern surveys, cosmic 
variance accounts for a large fraction of the uncertainty 
budget. Thus, we wish to quantify and compare the ef- 
fects of cosmic variance on the uncertainties in our mea- 
surements and those from the literature. First, we con- 
sider uncertainties due to cosmic variance in all recent 
measurements of the LF normalization in the literature. 

5.1. Cosmic Variance in Recent Measurements of the 
NIR LF normalization 

iDriver fc RobothamI (|2010[ ) have quantified cosmic 
variance for M* galaxies in the SDSS volume. They 
derive the following empirical formula to calculate the 
expected systematics due to cosmic variance in any low- 
redshift survey 



Cosmic Variance(%) = [1.00 - 0.03V(A/B) - 1] 
X (219.7 - 52.41ogio[AB x 291.0] 

-1-3.21 (Iogio[AB X 291.0])2)/v/NC/291.0(20) 

In the above formula, A and B are the transverse co- 
moving lengths at the median redshift of the survey, and 
thus provide an aspect ratio A/B for rectangular sur- 
veys (A > B) . C is the comoving length of the redshift 
interval surveyed, and N is the number of independent 
sightlines composing the survey. Given this formula, it 
is clear that surveys with multiple sightlines suffer less 
from cosmic variance systematics (by a factor of I/^/N) 
than single field surveys, and surveys with large aspect 
ratios are superior to those of square fields. 

We calculated the expected systematics due to cos- 
mic variance for each survey from the literature, given 
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depth, area, aspect ratio, ai id number of fields as input s 
to the empirical formula of [Driver fc RobothamI ()2010D . 
We were not able to account for the fact that, in many 
cases, these surveys consisted of several sightlines, each 
having a different size and aspect ratio, so these estimates 
are not meant to be extremely rigorous, but rather to 
provide a means of comparing the relative effects of cos- 
mic variance from one survey to the next. The results of 
this comparison are listed in Table [S] Below, we briefly 
discuss this estimate of cosmic variance in each survey 
from the literature. 

The first efforts using large sur veys to measure th e lo- 
cal NIR LF normalization were ICole et all ((200l and 
iKochanek et al.l ()2001l ). Both of t hese survey s used 
2MASS for their NIR photometry, and lCole etall (|2001[ ) 
used 2dFGRS spectroscopy for 5, 683 galaxies at a me- 
dian redshift z - 0.06, while IKochanek et al.l ()2001l ) 
used the ZCAT data f rom the CfA2 reds hift survey 
(IGeller fc Huchral IT989I: iHuchra et all [TOOl) for 3, 878 
galaxies at a median redshift z ~ 0.023. We find that 
both these surveys should suffer from systematics due to 
cosmic variance at the ~ 10% level. iCole et all ()2001l ) 
estimate cosmic variance systematics at the 15% level, 
using galaxy counts from 2MASS, and propagate this es- 
ti mate through the ir error estimation. 

lEke et all (I2005D provide an updated version of the 
study done bv 'Co le et all (|200lD using 2MASS and the 
completed 2dFGRS to include 15, 664 galaxies in the de- 
termination of the band LFs. Given the larger survey 
size, the estimated systematics due to cosmic variance in 
this sample are at the ~ 6% level. 

The la rgest study t o date of the local NIR LF was 
done by iJones et al.l H2006t ). They combine 2MASS 
photometry with the 6-deKre e Fiel d Galaxy Redshift 
Survey (6dFGS, Jon es et al.l 120041 ) to measure the 
J, if, and if— band LFs using ^ 60,000 galaxies at a 
median redshift z ^ 0.054. While this study includes 
more galaxies than previous works, the fact that it is rel- 
atively shallow implies that cosmic variance still persists 
at the 5% level. 

At med i an re dshifts z ~ 0.08, iBell et al.l (|2003l ) and 
iHiU et ahl (120101) proy i de a m easurement of the LF nor- 
malization. iBdl et al.l ()2003D use a sample 6, 282 galaxies 
selected from 2MASS with spectroscopy from the SDSS 
early data rele ase. They find a normalization quite sim- 
ilar to that of lEke et all (|2005f ). though with lar ger er- 
ror bars (cosmic varianc e at the ^ 10% level) due to 
the smaller sample size. iHill et all (|2010( ) combine the 
UKIDSS LAS, the SDSS, and t he Millennium Galaxy 
Catalog (MGC, iLiske et al.ll2003D to form a highly uni- 
form and 100% spectroscopically complete sample of 
~ 1800 galaxies in the J, ii, and if —bands. Their LF 
normalization agrees well with that of IBell et al.l (|2003[ ). 
though the relatively small volume of their sample im- 
plies possible systematics on the ^ 30% level due to cos- 
mic variance. 

iSmith et al.l ()2009() use a sample of ^ 40, 000 galaxies 
at a median redshift z ~ 0.1, drawn from the UKIDSS 
LAS with spectroscopy from the SDSS, to study the 
if-band LF. Their sample includes ah UKIDSS LAS 
galaxies at ifvoga < 16 that have secure redshifts in 
the SDSS. While this is not a spectroscopically complete 
magnitude-limited sample, it is one of the largest samples 
used to study the if —band LF and samples a large vol- 



TABLE 4 

NIR Luminosity Function Normalizations'' 

IN EACH OF six'' FIELDS FOR THE REDSHIFT 
RANGE 0.1 < 2: < 0.3 



Field 


X 10=^ 


0^ X 10^ 


X 10^ 


All Fields 


15.9±L7 


15.3±1.6 


14.6±1.7 


A370 


24.4±4.9 


23.1±4.7 


23.0±4.5 


CLANS 


14.4±3.5 


13.9±3.4 


12.7±3.1 


CLASXS 


12.0±3.3 


11.5±3.2 


10.3±2.8 


CDF 


10.9±6.4 


12.5±7.4 


10.4±6.3 


SSA13 


17.3±7.7 


15.5±7.7 


15.3±7.6 


SSA17 


13.3±5.6 


15.0±6.2 


N/A 



^ all 0* values arc given in units of IC" Mpc~^. 
These values were calculated using the 1 /Vmax 
method and fitting a Schcchter function with fixed 
a — 0.91 from applying STY in the — band, and 
fixed M* values (given in Table [3] derived in each 
bandpass using STY with a fixed a — 0.91. This 
method is described in Scction l4?4l 
^ only five fields in the case of the i^— band where 
the SSA17 field lacks photometry 

ume. For this survey size, we estimate cosmic variance to 
be at the 2%. However, given the relative spectroscopic 
incompleteness, other biases may be a significant source 
of systematic uncertainties. 

At median redshifts in the range 0.1 < z < 0.3, very 
few studies of the NIR LF have been done. Furthermore, 
the relatively small sample sizes imply that measure- 
ments of the LF normalization over this range are sub- 
ject to large systematic un certainties due to cosmic vari- 
ance. iHuang et al.l ()2003D use a sample of 1,056 galax- 
ies to study the if— band LF at a median redshift of 
z ~ 0.14. Their measured normalization is essentially in 
agreement with all measurements at lower redshift, given 
their statistical erro r s plus cosmic variance at the ~ 10% 
level. iFeulner et al.l (|2003[ ) use a sample of 210 galaxies 
to study the J and if —band LF at a median redshift of 
z ^ 0.2. The size of their survey implies uncertainties 
due to cosmic variance at the ~ 30% level, putting their 
measured normalization in agreement with all previous 
studies. 

Our survey samples the same redshift range as 
IFeulner et al.l (|2003( ) , but our larger area and deeper com- 
pleteness limits mean we have a factor of four more galax- 
ies in our study, and the expected uncertainties due to 
cosmic variance have been reduced by a factor of two. 
As noted above, our measured normalization is in good 
agreement with studies at z 0.1, but would appear to 
be in slight tension with studies at z ~ 0.05. In the fol- 
lowing section we explore the field to field variability in 
the LF normalization for our survey, and compare with 
simulations using the COSMOS Cone mock catalogs. 

5.2. Cosmic Variance in Our Survey 

We are measuring the normalization of the NIR LF 
over the redshift range 0.1 < z < 0.3, or comoving dis- 
tances of 300 - 850 h-^ Mpc. Thus, with a total of 
2 deg^ on the sky we are sampling the normalization over 
a volume of ~ 10^ Mpc'^. With typical large-scale 
structures existing on at least ^ 100 Mpc scales, we 
expect (and indeed find) significant field to field variation 
in the measured normalization. 

In Figure El we show the J, H, and if —band LFs for 
0.1<z<0.3ona field to field basis. In all cases the 
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Fig. 8.— The J, H, and ii'-band LFs of galaxies for 0.1 < 2 < 0.3 
in the J(a), -ff(b), and K{c) bands. The solid black line shows 
Schechter fits to the combination of all six fields (black asterisks, 
only five fields in the ii"— band where SSA17 lacks photometry). 
The colored lines show Schechter fits for each of the six fields in- 
dividually. These LFs are derived in the same way as those pre- 
sented previously, with Schechter fits derived from the results of 
the l/Vmax method using a fixed a and M* derived from t he STY 
method. Error bars show the 1 a Poisson errors from TJeEreli 
l|1986). The Schechter function normalizations for each individual 
field in this redshift selection are given in Table U] We note that 
the normalization in the A370 field is roughly a factor of two higher 
than in the other fields, even though the A370 cluster itself, at a 
redshift of z = 0.37 is excluded. 
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Fig. 9. — (a) Redshift versus ii"— band absolute magnitude in 
an example selection of galaxies {K < 18) for the combination of 
six sightlines of 1/3 deg'^ each taken at random from the COSMOS 
cone mock catalogs, (b) A histogram of all redshifts in the six mock 
catalogs. The solid black line shows the total over six fields and 
the colored lines represent the six fields individually. We note the 
similarity in the distribution of redshift and absolute magnitudes 
in these example mock catalogs in comparison to our observations 
shown in Figure [Tl 



A370 field shows a much higher normalization than the 
other five fields, even though the A370 cluster itself, at 
redshift z ~ 0.37, is not included in these LFs. The solid 
black line shows Schechter fits to the combination of all 
six fields. The colored lines show Schechter fits for each 
of the six fields individually. 

In Table S] we give the normalizations for each field 
and the combination of fields in the range 0.1 < z < 0.3. 
These values were calculated using the 1/Vmax method 
and fitting a Schechter function with fixed a = 0.91 
from applying STY in the if— band, and fixed M* val- 
ues (given in Table [3] derived in each bandpass using STY 
with a fixed a = 0.91). 

We find agreement between the normalizations in all 
fields except the A370 field, where the normalization is a 
factor of '-^ 2 higher. Upon investigation of this apparent 
excess, we find that two sheet-like overdensities exist in 
the A370 field at redshifts of z - 0.18 and z ~ 0.25 (both 
clearly visible in Figure [IJ. These overdensities are very 
narrow in redshift space, spanning the redshift ranges 
0.17 < z < 0.19 and 0.245 < z < 0.255, respectively. 
To determine whether our survey appears to be consis- 
tent with cosmological simulations, and to further ex- 
plore the expected systematics due to cosmic variance in 
this study, we use the COSMOS Cone mock catalogs. 

5.3. Comparison With The COSMOS Cones 

The COSMOS cones are composed of 24 pencil 
bea m mock ca talogs crea ted by Manfred Kitzbich- 
ler (|Kitzbichler k White .20071) for the Cosmic Evo- 
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Fig. 10. — Percent deviation from the mean NIR LF normal- 
ization for simulated surveys of 6 (black), 12 (red), and 24 (blue) 
sightlines. The histograms show the distribution in measured nor- 
malization for 1000 realizations of the simulated surveys from the 
COSMOS cone mock catalogs. The dashed lines show a Gaus- 
sian fit to each histogram. Each survey is intended to mimic our 
observed data, in that the simulated data are magnitude limited 
(K < 18) and restricted to the redshift range 0.1 < ^ < 0.3. 
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Fig. 11. — Systematic uncertainty in percent (1 c) due to cosmic 
variance as a function of the number of sightlines in surveys com- 
posed of multiple square sightlines of 1/3 deg^ each {K < 18 and 
O.l < z < 0.3. The solid line shows the empirical estimate from the 
formula of Driver & Robotham (2010) and the dashed line shows 
our estimate calculated from the COSMOS Cone mock catalogs. 



lution Survey collaboration (COSMOS, iScoville et all 
l2007bD and made public in the Millennium Database 

onlin43. These catalogs were created from the 

Millenium simulation (Springel et al. 2 0051) us i ng the 
galaxy formation algorithm of iDe Lucia et al.l (I2OO60 
and th e merger tree prescriptions of lOe Lucia fc BlaizotI 
(j2007[) . The COSMOS cone catalogs have been used 
primarily to study galaxy formation and evolution, 
as well as large-scale structure, at relatively high 
redshifts iS covillc et aT 2007a|; 
iKnobel et a l. 2009; Kovac etaLl 



McCracken eFaLl 120071: 
2010D. We use the cat- 



alogs here to compare with our observations of the NIR 
LFs at 0.1 < z < 0.3. 

The COSMOS cone mock catalogs represent 8 sight- 
lines from each of 3 different origins within the Millenium 
simulation. Each of the COSMOS cone pencil beams 
covers an area of 1.4 x 1.4 degrees, or ~ 2 deg^. Thus, 
each mock catalog is roughly the equivalent area of our 
entire survey, and so we divide each mock catalog into 
six separate sightlines. This gives us a total of 144 mock 
sightlines to compare with our observations. /iT— band 
apparent magnitudes are provided in the mock catalogs 
(though not J or band), so we take a magnitude lim- 
ited sample from the mocks of K < 18 to compare with 
our data. 

From the 144 sightlines, we first select six at ran- 
dom. We show an example of the distribution in ab- 
solute iiT— band magnitude and redshift for six random 
sightlines selected from the COSMOS cones in Figure |9l 
This apparent magnitude selection of i^T < 18 typically 
generated a catalog of ^ 1500 objects in the combination 
of six sightlines. We note the similarity in the distribu- 
tion of redshift and absolute magnitude in comparison to 
our observed data for iJ < 18 shown in Figure [TJ 

Next, we construct LFs to study the variability in 
the normalization from field to field. Rather than us- 
ing the absolute magnitudes provided in the catalogs, 
we attempted to mimic our methodology applied to the 
observed data by using the apparent A'— band magni- 
tudes from the mock catalogs, and converting to abso- 
lute magnitude via a distan ce modulus, a ii'— correction 
from lMannucci et al.l ()2001[ ). and an E{z) = z correction. 
We then compared the normalization of the LFs in six 
randomly selected mock catalogs in the same manner as 
with the observational data, where we estimated a and 
M* with the STY method, then found the normalization 
fitting a Schechter function to the l/Vmax results with a 
fixed a and M*. 

We completed 1000 realizations of the process of select- 
ing six sightlines at random and evaluating their LFs. We 
found that we have a ~ 25% chance of drawing a survey 
like ours, with one overdense sightline, from the mock 
catalogs. Thus, while the appearance of overdensities 
like those in the A3 70 field might not be characteristic of 
sightlines chosen at random from simulations, they are 
not strongly atypical. However, we note that in both ob- 
servations and simulations, massive galaxy clusters tend 
to lie at the confluence of large filamentary structures, 
such that choosing to observe in the direction of a galaxy 
cluster may increase the chances of finding overdensities 
along the line of sight. 

Next, we measured the expected systematic uncer- 
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tainty due to cosmic variance by looking at the distri- 
bution of measured normalizations for our 1,000 mock 
surveys of six sightlines. We find the distribution is well 
fit by a Gaussian (see Figure [T0|) with a 1 a variability 
of ^ 13% in the measured normalization over the 1, 000 
combinations of six random sightlines. 

We repeated the process of 1000 realizations of sim- 
ulated surveys with the number of sighthnes (1/3 deg^ 
each) ranging from 1 — 30. Also shown in Figure [10] is the 
distribution in measured normalization for surveys of 12 
and 24 sightlines for comparison. 

In Figure [Tl] we show the 1 a systematic uncertainty 
(in percent) due to cosmic variance as a function of the 
number of square 1/3 deg^ sightlines. The dashed line 
shows our results from simulated surveys drawn from the 
COSMOS cones. The expected systematics gi ven the em- 
pirical formula of lDriver fc RobothamI ()2010D are shown 
as a solid line. We find that the two estimates agree 
rather well, and hence we confirm their results derived 
from observed data with our analysis of the COSMOS 
cones. Given these results, we find that doubling the 
number of sightlines in our survey should reduce the ef- 
fects of cosmic variance to the ~ 10% level, while qua- 
drupling the survey should reduce these systematics to 
the ~ 5% level. 

6. SUMMARY 

An accurate measurement of the average space den- 
sity of galaxies just beyond the local universe (z > 0.1) 
is of interest, in particular, because recent cosmological 
modeling efforts have shown that if the matter density 
at z < 0.1 is low by roughly 50%, then the apparent ac- 
celeration of the expansion of the universe observed in 
type la supernovae could simply be a p roduct of our lo- 
catiqn with respect to local st ructure (jAlexander et al.l 
[200l IBoleiko fc Sussma J'20TT). 

We have presented a study of the NIR LF normaliza- 
tion at redshifts of 0.1 < z < 0.3. This study is based 
on a highly complete (> 90%) spectroscopic sample of 
812 galaxies, selected to be brighter than 18*'' magni- 
tude in the band, over six widely-separated fields at 
high Galactic latitudes. While a relatively small sample 
of galaxies, our survey represents a factor of four increase 
in the number of galaxies in previous studies of the NIR 
LF in this redshift range, which implies a factor of two 
decrease in the systematic uncertainties due to cosmic 
variance. 

We construct the NIR LFs for each field individually 
and for the combination of all six fields. We compare 
the product of the normalization and characteristic lu- 
minosity ((j)* L*) in the combination of our six fields with 
lower redshift measurements to test for a higher lumi- 
nosity density at z > 0.1. We find our measurement of 
(j3*L* to be in agreement with other studies at median 
redshifts near z ^ 0.1, but roughly 30% higher than the 
error-weighted mean at (z) 0.05. Our measurement of 
(j)*L* for 0.1 < z < 0.3 could be considered a conserva- 
tive underestimate of the true value of cf)* L* because we 
avoid known galaxy clusters in this redshift range. 

We use the CO SMOS cone mock catalogs, a nd the em- 
pirical formula of lDriver fc RobothamI ()2010D to investi- 
gate the effects of cosmic variance on our measurement 
of the NIR LF normalization and on those from the liter- 
ature. Despite the uncertainties due to cosmic variance. 



it would appear that data from the literature are begin- 
ning to suggest a rising luminosity density from z = 0.05 
to z = 0.1. While this conclusion is tentative at best, it 
would not be terribly surprising, as a rising LF normal- 
ization over this redshift range would be predicted given 
the results of several NIR galaxy counts surveys over the 
past 10 — 20 years. 

More importantly, however, we require knowledge of 
the luminosity density at z > 0.1 to understand whether 
or not the local universe is under-dense. We find that 
measurements of the NIR LF at z > 0.1, including our 
own, are still too uncertain to provide a robust compar- 
ison with lower redshift measurements. While inconclu- 
sive, this result is noteworthy given the possible implica- 
tions for locally measured cosmological observables. 

Surveys that are currently underway, including our 
own campaign of spectroscopic follow-up on fields that 
have existing NIR photometry, will soon provide a more 
robust answer to the question of whether or not we reside 
in a large local underdensity. 
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